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Longitudinal optical phonons with oxygen character were measured in Lao.rSro.sMnOs by inelas- 
tic neutron scattering in the (10 0) cubic direction and results were compared with shell model 
predictions. Measurements were performed in several Brillouin zones, which enabled us to identify 
the eigenvectors independent of the shell model. All major disagreements between model predic- 
tions and experimental results are primarily due to the anomalous downward dispersion of the 
bond-stretching vibration. The main new result is that the rhombohedral distortion of the cubic 
lattice makes the bond-stretching vibrations interact strongly with bond-bending modes folded into 
the cubic Brillouin zone. 



I. INTRODUCTION 

Interplay between the charge, spin, and lattice 
degrees of freedom in perovskite manganites re- 
sults in a multitude of unconventional properties 
of fundamental as well as practical interest. Fer- 
romagnetic alignment of Mn core spins appears at 
some stoichiometrics due to the double exchange 
interaction and the ferromagnetic-paramagnetic 
(FM) transition in these systems is accompanied 
by large magnetoresistance (MR) . As predicted by 
theory^, polarons forming and condensing above 
the FM transition temperature, Tc, have been ob- 
served by neutron scattering in Lao.TCao.sMnOa 
and other manganites exhibiting colossal magne- 
toresistance (CMR )^-^-i .. Trapping of conduction 
electrons by these polarons is held responsible for 
the CMR effect. While the presence of polarons 
above Tc demonstrates that electron-phonon cou- 
pling due to the Jahn- Teller effect is crucial to un- 
derstanding the essential physics, relatively little 
is known about phonons in these materials. 

"113" manganites such as Lai_a;Caa;Mn03 and 
Lai-a^Sra^MnOa have a cubic perovskite structure 
except for a rotation of the MnOe octahedra, which 
doubles the unit cell.^ In this paper we will use 
the notation based on the cubic unit cell of the 
undistorted structure, which makes the analysis 
easier to understand. Use of this notation re- 
quires folding in extra phonon branches that arise 
from the doubled unit cell. (For example, some 
phonons observed only at the zone boundary in 



the undistorted structure have nonzero structure 
factor at the zone center in the distorted struc- 
ture, etc.) A previous study of oxygen phonons in 
LacrSrcsMnOa^-^ found anomalous downward dis- 
persion of the bond-stretching vibration in the (1 
0) direction, but it did not examine effects of 
interactions of these vibrations with the folded-in 
branches allowed by the tilt of the octahedra. 



According to the presently accepted view of 
the physics of the ferromagnetic manganites, the 
strongest temperature dependence is expected of 
the Jahn- Teller modes, whose wavevectors are (0.5 
0.5 0.5) and (0.5 0.5 0). But investigations of 
LacrCacsMnOs^ and Lao.rSro.sMnOs^ found that 
these phonons did not have strong temperature 
dependence. Instead, a significant intensity loss 
on raising the temperature towards the FM tran- 
sition was reported for phonons with wavevec- 
tors dispersing in the 10 direction. These re- 
sults motivated us to go beyond the study of Re- 
ichardt and Braden.^ We measured phonon dis- 
persions with higher resolution where possible and 
identified eigenvectors of q=(X 0) phonons in 
Lao.TSro.sMnOa including interactions with the 
folded-in branches into the analysis. Although the 
crystals are twinned, the (10 0) direction is the 
same for the two domains. 
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II. EXPERIMENTAL AND 
CALCULATION DETAILS 

The experiments were carried out on the triple- 
axis spectrometer IT located at the ORPHEE 
reactor using doubly focusing Culll and Cu220 
monochromator crystals and PG002 analyzer fixed 
at 13.7, 14.8, or 30.5 meV. Our samples were two 
high quality single crystals of Lao.rSro.sMnOs with 
the FM transition temperature measured at 355K. 
The volume of each crystal was O.Scm^. The re- 
sults were the same for the two crystals. All mea- 
surements were performed in the ferromagnetic 
phase at lOOK (we did not go to lower tempera- 
tures for technical reasons). Since the FM transi- 
tion temperature is much higher, we expect to get 
the same results at lower temperatures. This was 
confirmed by repeating some scans at 12K. 

We performed measurements in several Brillouin 
zones and, guided by shell model predictions and 
symmetry constraints, found by trial and error ap- 
proximate phonon eigenvectors that correctly give 
the observed intensities. In this procedure we ig- 
nored small transverse atomic displacements, since 
these would not affect longitudinal phonon intensi- 
ties in the Brillouin zones that we measured. They 
should exist due to the rhombohedral distortion. 



III. CALCULATION RESULTS 

Ignoring the rhombohedral distortion of the lat- 
tice, we expect to observe two optical oxygen 
branches intrinsic to the cubic perovskite struc- 
ture above 35 meV. These are the bond-stretching 
and bond-bending modes. The rhombohedral dis- 
tortion doubles the unit cell along one of the (111) 
directions and thus the number of phonon branches 
is expected to double as well. Shell model calcu- 
lations show that the bond-bending and the bond- 
stretching branches should survive in the distorted 
structure, plus two additional branches of bond- 
bending character should be folded into the cubic 
Brillouin zone due to the unit cell doubling. At 
the zone center one of the folded-in modes has 
a vanishing structure factor, whereas the other 
one has a vanishing structure factor at the zone 
boundary. Thus there should be three modes with 
nonzero intensity at the zone center and the zone 
boundary and four modes in the middle of the 
zone. According to the model, the original bond- 
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FIG. 1: Longitudinal oxygen phonon dispersions 
along the 10 cubic direction calculated using the 
shell model. Only the branches with nonzero structure 
factors are shown. 



bending branch should disperse upwards, whereas 
the bond-stretching and the folded-in branches 
should be approximately flat. (See Fig. 1) 

Structure factors of the bond-bending and 
folded-in modes are very sensitive to the value of 
Mn-O-Mn bond angle whereas the structure factor 
of the bond-stretching vibration is not. The bond 
angle depends on the magnitude of the rhombo- 
hedral distortion of the lattice and is reported to 
be 164° for Lao.rSro.sMnOa.^ In the absence of the 
rhombohedral distortion, i.e. when the Mn-O-Mn 
bond angle is 180°, all longitudinal phonon inten- 
sities should be proportional to Q^. 



IV. EXPERIMENTAL RESULTS 

A. Zone Center 

We observe the zone center phonons at 42meV, 
53meV, and 72meV. (Fig. 2a) According to the 
shell model these should be the energies of the 
bond-bending, folded-in, and the bond-stretching 
modes respectively (Fig. 2b). This assignment can 
be verified by structure factor calculations for the 
corresponding polarization patterns. The 42 meV 
bond-bending vibration decreases in intensity by a 
factor of 1.82 (Fig. 2a) while the structure factor 
of the corresponding (leftmost) polarization pat- 
tern in figure 2b gives 2.17. The increase of the 
53meV phonon intensity from Q=(4 0) to Q=(5 
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FIG. 2: (Color online) a) Zone center phonons mea- 
sured in two Brillouin zones. Resolution was lower at 
Q=(5 0). b) Approximate eigenvectors of the zone 
boundary modes at lOOK deduced from the measured 
intensities. Mode frequencies left to right are: 42mev, 
53meV, 73meV. Only pairs of MnOe octahedra are 
shown. 



0) is 1.42 vs the calculated value of 1.96. The 
72meV mode intensity increases by a factor of 1.28 
vs. the calculated 1.53. Considering that we ignore 
small deviations from the cubic high symmetry di- 
rections as discussed above as well as small admix- 
tures of other modes, structure factor calculations 
confirm our assignment. 
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FIG. 3: (Color online) a) Zone boundary phonons 
measured in four Brillouin zones, b) Approximate 
eigenvectors of the zone boundary modes at lOOK in 
the MnOe octahedra deduced from the measured in- 
tensities. Mode frequencies left to right are: 38mev, 
47meV, 61meV. 



B. Zone Boundary 

q=(0.5 0) (zone boundary) is another high 
symmetry point of the reciprocal lattice. We per- 
formed measurements in different Brillouin zones 
in order to gain insight into the zone boundary 
phonon eigenvectors. Three longitudinal oxygen 



modes appear at 38, 47, and 61 meV. 

We begin the analysis with the 61 meV phonon. 
Its intensity decreases with increasing Q and al- 
most vanishes at Q=(5.5 0). That is pre- 
cisely what is expected from the bond-bending 
zone boundary mode (rightmost polarization pat- 
tern in figure 3b) as discussed for the zone center 
modes. In fact structure factor calculations for this 
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FIG. 4: (Color online) Comparison of the measured 
and calculated phonon intensities at the zone bound- 
ary for the 38meV, 47meV, and 61meV phonons. The 
calculated values are the predicted intensities assum- 
ing the corresponding polarization patterns in Fig. 3b. 
An arbitrary scaling factor (the same for all points in 
the three graphs) was applied to the calculated values. 



phonon eigenvector are in good quantitative agree- 
ment with the observed intensity at 61meV (Fig. 
4). The shell model predicts the correct frequency 
for this mode, thus we can assign it with a very 
large degree of confidence. 

The 47 meV phonon is the strongest at Q=(4.5 
0); it is weaker at Q=(2.5 0) and is absent at 
Q=(3.5 0) and Q=(5.5 0) (Fig. 3a). Structure 
factor calculations based on the middle polariza- 
tion pattern in figure 3b correctly reproduce the 
observed intensity at all measured Q. (Fig. 4) The 
47 meV phonon is mostly bond-stretching in char- 
acter, but with a significant bond-bending com- 
ponent of folded-in character (The bond-bending 
component in the 47 meV mode is orthogonal to 
the bond-bending vibration at 61 meV). 

The 38 meV phonon appears most strongly at 
Q=(3.5 0), it is weaker at Q=(5.5 0) and 
Q=(4.5 0), and is very weak at Q=(2.5 0) (Fig. 
3a). The leftmost polarization pattern in figure 
3b gives the correct intensities at all measured Q. 
(Fig. 4) Just like for the 47 meV phonon, its eigen- 
vector is a superposition of the bond-stretching 
and folded-in bond-bending vibrations, but here 
the bond-bending vibrations have a larger ampli- 
tude. 

The coupling between the bond-bending com- 
ponent and the bond-stretching one in the 38 and 
47meV modes is exclusively due to the tilt of the 
octahedra and should decrease with increasing Mn- 
0-Mn bond angle. Most of the bond-stretching 
character at the zone boundary is at 47 meV with a 



smaller contribution at 38 meV. These energies are 
much lower than the prediction of the shell model 
and indicate strong electron-phonon coupling. 



C. Between the Zone Center and the Zone 
Boundary 

Deducing the eigenvectors between the zone cen- 
ter and the zone boundary as was done at the zone 
center and the zone boundary is much more diffi- 
cult due to lower symmetry. However, on a qual- 
itative level, this analysis is still useful. Figure 5 
shows the phonon spectra at three reduced recipro- 
cal lattice vectors each measured in two Brillouin 
zones. There are four phonon peaks in agreement 
with the predictions of the shell model. 

We start with the bond-bending branch, which is 
the easiest one to identify by its intensity decrease 
towards large Q (see above). The intensity of the 
highest energy mode at X=0.4 and 0.3 decreases 
from Q=(4-X 0) to Q=(4+X 0). (Fig. 5) It is 
also very strong at Q=(2.3 0) (data not shown). 
Such behavior is similar to the 61 meV mode at 
the zone boundary, thus we conclude that the high- 
est energy phonons at X=0.4 and 0.3 are of bond- 
bending character. The bond-bending branch orig- 
inates at 42 meV at the zone center (Fig. 2) and 
disperses upwards crossing the 53 meV branch at 
q^{0.15 0)(data not shown). Mixing with the 
bond-bending branch may be the reason that the 
53meV peak is stronger at Q=(3.8 0) than at 
Q=(4.2 0). Between q=0.1 and q=0.3 it hy- 
bridizes with the other branches emerging again 
as a pure mode at q=(0.3 0) - (0.5 0) around 
60 meV. 

Identification of the dispersion of the bond- 
bending branch allows us to assign the other 
modes to either the bond-stretching or the folded- 
in branches. The mode at q=(0.2 0) and 65 meV 
has mostly bond-stretching character, since its in- 
tensity is stronger at Q=(4.2 0) than at Q=(3.8 
0). We also know that the mode at q=(0.3 0) 
and 57meV is bond-bending. This indicates that 
the bond-stretching branch has a very steep down- 
ward dispersion between q=(0.2 0) and (0.3 
0). This steep dispersion is responsible for the 
very large energy width of this mode at q=(0.2 
0) due to the finite q-resolution of the spec- 
trometer. At q= (0.3-0.5 0) the bond-stretching 
mode continues to disperse to lower energies, but 
not as steeply. In this range of wavevectors it hy- 
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FIG. 5: (Color online) Phonons between the zone cen- 
ter and the zone boundary measured in two Brillouin 
zones 



bridizes with the two folded-in branches. As a re- 
sult we observe three peaks below 55 meV. All are 
of mixed folded-in/bond-stretching character as is 
clear from strong intensity changes between Q=(4- 
X 0) and Q=(4+X 0) (see discussion above for 
the zone boundary). Figure 6 illustrates the dis- 
persion of the bond-stretching vibration (marked 
with arrows) in the data covering Q= (4-4.5 0). 
It disperses downward away from the zone center 
and at q>(0.2 0) it contributes to more than one 
mode due to its interaction with other branches. 



V. SUMMARY 

Figure 7 shows the dispersion of the longitudinal 
oxygen phonons in the 10 direction summariz- 
ing the above results. The shell model gives a good 
prediction for the frequencies of the bond-bending 
and folded-in branches, however, it does not pre- 
dict a strong downward dispersion of the bond- 
stretching branch away from the zone center as pre- 
viously reported.^ This anomaly has been observed 
in almost all metallic perovskite oxides and is gen- 
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FIG. 6: Evolution of the phonon spectra between 
Q=4 and Q=4.5 0. Arrows mark modes with 
bond-stretching character. 




FIG. 7: (Color online) Circles represent measured lon- 
gitudinal phonon peak positions at lOOK. Vertical bars 
represent peak widths. Solid lines are guides to the eye 
showing dispersions of the bond-stretching (downward- 
dispersing) and bond-bending (upward-dispersing) vi- 
brations in the absence of interactions with other 
branches and each other. Their crossing of each other 
and of the folded-in branches (see text) results in com- 
plex phonon polarization patterns. Dashed lines show 
results of the shell model calculations. 



erally interpreted as a signature of strong electron- 
phonon coupling. As a result of the downward 
dispersion, the bond-stretching branch crosses the 
bond-bending branch at q=0.2 and hybridizes with 
the folded-in branches at q=(0.2-0.5 0). The 
bond-bending branch has an upward dispersion in 
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good agreement with the shell model. The folded- 
in branches are nearly fiat, but they are affected 
by the interaction with the bond-stretching branch 
at q=(0.2-0.5 0). 

VI. CONCLUSIONS 

In conclusion, we have demonstrated that the 
rhombohedral distortion of the lattice has a pro- 
found effect on the longitudinal oxygen phonons 
dispersing in the 10 cubic direction. It 
results in two nearly fiat folded-in branches 
which couple strongly to the upward-dispersing 



bond-bending and the downward-dispersing bond- 
stretching branches. The latter appears as a pure 
mode only near the zone center hybridizing with 
the folded-in branches throughout most of the zone 
including the zone boundary. 
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